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We consider metric-affine scenarios where a modified gravitational action is sourced by electrovac¬ 
uum fields in a three dimensional space-time. Such scenarios are supported by the physics of crys¬ 
talline structures with microscopic defects and, in particular, those that can be effectively treated 
as bi-dimensional (like graphene). We first study the case of f(R) theories, finding deviations near 
the center as compared to the solutions of General Relativity. We then consider Born-Infeld gravity, 
which has raised a lot of interest in the last few years regarding its applications in astrophysics and 
cosmology, and show that new features always arise at a finite distance from the center. Several 
properties of the resulting space-times, in particular in presence of a cosmological constant term, 
are discussed. 
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I. INTRODUCTION 

One hundred years after the development of Einstein’s 
theory of General Relativity (GRb many experimental 
results have confirmed its success Jl] . This theory, how¬ 
ever, is still troubled by the presence of space-time sin¬ 
gularities 0, where it breaks down and loses its reliabil¬ 
ity. In this sense, there is the widespread belief that the 
solution to this problem will be achieved once GR and 
quantum theory merge into a single, unified framework 
-a quantum theory of gravity-, like in string theory [3] or 
loop quantum gravity. Another approach to this problem 
consists of obtaining some insights using classical theo¬ 
ries of gravity including higher-order corrections in the 
curvature invariants [5j] • The simplest such example is 
that of f(R) theories, where / is a given function of the 
curvature scalar, R. Indeed, such an approach has al¬ 
lowed to shape different astrophysical and cosmological 
phenomena (see e.g. [5] for a review). 

Given the mathematical complexity of GR and its ex¬ 
tensions, a useful way to study their physics relies on con¬ 
sidering their three dimensional counterparts Q. Here 
analytic solutions can be obtained and characterized, 
which might provide useful insights for four dimensional 
black hole physics. Indeed, three dimensional black holes 
may arise from effective descriptions of four dimensional 
gravity in string theory Q. In this simplified scenario, a 
relevant finding is that of the BTZ solution Q (see also 
I) , where a nonsingular Anti-de Sitter space-time arises 
as a bound state separated from the continuous black 
hole spectrum by a mass gap. This solution and its exten¬ 
sions have been widely studied due to their potential im¬ 
plications within the AdS/CFT correspondence 0 and 
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systems of condensed matter physics. Indeed, regular 
solutions sourced by an electromagnetic field date back 
to the Bardeen solution 0 and its nonlinear extensions 
[jjj. In three space-time dimensions, regular black hole 
solutions have been obtained as generalizations of the 
Bardeen solution [13]. Other three dimensional space- 
times include wormhole solutions in nonlinear electrody¬ 
namics mm , charged shells scenarios Ifil. b lack holes 
with dilaton fields 17j, magnetic solutions [lq], and fur¬ 
ther extensions of the BTZ solution 0. 

When dealing with extensions of GR, however, one is 
immediately faced with the troubles of higher-order field 
equations and ghost-like instabilities. One is then forced 
to constrain the parameters of the theory (when possible) 
in order to get rid of such undesirable features (20]. How¬ 
ever, another approach to this problem comes from the 
study of crystalline structures and, in particular, Bra- 
vais crystals in solid state physics (2lj| . At the micro¬ 
scopic level they are defined by a discrete net-like ar¬ 
rangement of atoms. However, it has been found that at 
the macroscopic level they admit a description in terms of 
an effective differential geometry [22j . If the crystal con¬ 
tains no defects on its microscopical structure, then the 
appropriate geometry is the standard Riemannian one 
of GR. However, all crystals contain defects of different 
sorts, and the Riemannian description no longer holds. 
A metric-affine geometry becomes necessary [23| . If such 
defects are point-like then one has non-metricity [24| (the 
line-like defects called dislocations require instead Car- 
tan’s torsion 0), which implies that metric and affine 
connection are not constrained by the standard compati¬ 
bility condition of GR (i.e., the latter is not given by the 
Christoffel symbols of the former), but are instead inde¬ 
pendent entities (see [26] for a pedagogical discussion). 
Independent variations of the action are then performed 
with respect to the metric and the connection, and the 
specific form of the latter is determined through the cor¬ 
responding equation • It turns out that this approach 
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leads to second-order equations and absence of dynamical 
instabilities for a large family of gravitational theories. 

The application of this framework to a three dimen¬ 
sional space-time could be of great relevance regarding 
some systems of solid state physics like graphene, which 
consists of an extremely thin layer of hexagonal atoms 
[28|, so it can be effectively treated as bi-dinrensional. 
The presence of defects on its nricrostructure raises the 
interest in the consideration of metric-affine theories. Let 
us point out, however, that these geometrical scenarios 
have received little attention so far and, consequently, 
they are still poorly understood. The goal of this paper 
is to made a modest step in this direction, and consider 
metric-affine gravities with non-metricity in a three di¬ 
mensional space-time. We focus on two classes of theo¬ 
ries, f(R) and Born-Infeld gravity, the latter being sup¬ 
ported by some results of solid state physics systems with 
defects. The dynamics of the new gravitational system 
is excited by considering an electromagnetic (Maxwell) 
field. 


II. f(R) THEORIES IN METRIC-AFFINE 
APPROACH 

The action of our theory is given by 

S= 2^l d ' x ^~3f{R) + SmtevuM, (!) 

where n 2 is Newton’s constant in suitable dimensions, g 
is the determinant of the space-time metric g^y, f(R) is a 
given function (unspecified at this stage) of the curvature 
scalar R = g llv R IJIJ with the definition R^ = R a ^, av , 
where 

R a P „ = + r“ A r^ - r“ A r^ (2) 

is the Riemann tensor. S m is the matter action and 
it depends on the matter fields, denoted collectively as 
ipm, and is assumed to only couple through the met¬ 
ric. In these definitions, the connection T*„ is a priori 
independent of the metric (metric-affine or Palatini ap¬ 
proach). In addition, we shall assume vanishing torsion, 
- T^ = 0 HH. As we will be interested 
in electrovacuum space-times, and to work as general as 
possible, the matter sector of our theory is that of non¬ 
linear electrodynamics, with action 

S m = j d 3 Xy/^gip(X), (3) 

where <p(X) is a given function of the field invariant 
X = —\F IIV F I1V constructed with the field strength ten¬ 
sor F^ u = d^Ay — d y A M , with A p the vector potential. 
To obtain the field equations for the theory defined by 
Eqs.© and ([3]) with the constraints above, we perform 


independent variations with respect to the metric and 
the connection, which yields (the details of this variation 
can be found in [30] for arbitrary dimension, and thus we 
only bring here the final result) 

IrR^u ~ ^9^ = k 2 T^ v , ( 4 ) 

v{{V=gfR9n = 0, (5) 

where we have introduced the short-hand notation /r = 
df/dR , and T pu = — ^=|^- is the energy-momentum 
tensor of the matter, which is obtained as 

T pv = 2ip x F pa F a y — ip(X)g fl y, ( 6 ) 

where <px = dtp/dX. Non-metricity appears explicitly 
in this framework through Eq.©. Before solving these 
equations, let us first take a trace in d4j), which yields 

RfR -\f = k 2 t, (7) 

where T = 2p — 4 Xtpx is the trace of T^y. Note that 
© is not a differential equation but just an algebraic 
relation R = R{T) extending the GR one, R = —k 2 T, 
to the case of metric-affine f(R) gravity coupled to non¬ 
linear electrodynamics. This implies that f(R) is just a 
function of the matter fields. 

Now, to solve the system of field equations © and © 
we take advantage of the fact that m can be formally 
rewritten as the standard metric-compatibility condition, 
namely 

= 0, (8) 

which means that is the Levi-Civita connection of 
a rank-two tensor /i M „. A little algebra shows that the 
relation between g P y and h pv is given by the conformal 
transformation 

V = ^ = /* V". ( 9 ) 

In terms of h^y the metric field equations © admit the 
following simple representation 

= + (10) 

This is a set of Einstein-like second-order field equations 
for the metric h^y , with the right-hand-side only depend¬ 
ing on the matter fields. The field equations for g^y fol¬ 
low from those of h pv through the matter-dependent con¬ 
formal transformation of Eq.© and, therefore, they are 
second-order as well. In vacuum, T fJ v = 0, the field equa¬ 
tions boil down to those of GR with possibly a cosmo¬ 
logical constant term, depending on the particular f(R) 
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function chosen. In summary, the system of equations 
(Hod, with the relations ©, 0 and 0 provide a solu¬ 
tion for a given problem once the forms of f(R) and <p(X) 
are specified. 

We note that the field equations 03 are in sharp con¬ 
trast with the more standard case of the metric formula¬ 
tion of f(R) gravity (see e.g. (31, 32|] for reviews of the 
theory and confrontation between these two approaches), 
where they take the form 

- ^g^ - V^V^/r + = k 2 T ^ (11) 

This is a set of fourth-order differential field equations, 
whose resolution is, in general, a non-trivial task. In this 
sense, simplifications such as that of constant curvature, 
R = Ro, which allows to remove the contribution of the 
fourth-order terms, allows for the obtention of explicit 
solutions (see e.g. [3l|). In our case, the presence of non- 
metricity, = V ^g v \, which measures the failure of 

the independent connection to be metric with respect to 
g^v [33j], introduces a relative (conformal) deformation 
between the Riemannian structure associated to and 
that associated to the physical metric g^ u . Note that, 
because of the different structure of the field equations 
and methodology followed to find analytical solutions in 
the metric and metric-affine approaches, here we will be 
mostly interested on comparing the deviations with re¬ 
spect to GR. 

A. Electrovacuum solutions 

We are interested in solving the field equations for a 
three dimensional, static, spherically symmetric line ele¬ 
ment of the form 

ds 2 = gtt.dt 2 + g rr dr 2 + r 2 dd 2 . (12) 

For a purely electric field the only-nonvanishing com¬ 
ponent of the field strength tensor in this setting is 
F tr ^ 0. From the matter field equations, SS m /5A v = 
0 -A V f i(y/—g<fixF IJ ' 1 ') = 0, it follows that this compo¬ 
nent satisfies 

F tr = , (13) 

fy/-gttgrrV-X 


To proceed further, we need to specify the choice for 
f(R) and <p(X). For the gravity part, we consider the 
theory 

f(R) = R + aR 3/2 (15) 

by computational convenience, since the trace equation 
0 yields the linear relation R = —2 k 2 T, thus simplify¬ 
ing the analytic resolution of the field equations. In this 
theory a is some constant (with dimensions), assumed 
to be small, which controls the deviation with respect 
to GR solutions, to which the new configurations will be 
compared. For the matter section, we take the choice 
<p(X) = X/(47t), corresponding to Maxwell electrody¬ 
namics, which in three dimensions yields a non-vanishing 
trace T = q 2 /( 47t?’ 2 ) (in four dimensions, however, the 
Maxwell energy-momentum tensor is traceless and one 
needs to use a nonlinear theory of electrodynamics like 
the Born-Infcld one in order to obtain modified dynamics 
[51). Thus we have R = 2 r 2 /r 2 , where r 2 = k 2 <? 2 /(47t) 
is a charge scale. Next we introduce two line elements, 
one for the metric h as 

ds 2 = —A(x)e 2 ^ x ' > dt 2 + ^ ^ dx 2 + r 2 (x)d0 2 , (16) 

and another for g^ as 

ds 2 = —B(x)dt 2 + C{x)dx 2 + r 2 {x)d6 2 , (17) 

where r 2 (x) and r 2 (x) are, in general, functions of the 
coordinate x, while £(x) and A(x) are functions to be de¬ 
termined by solving the field equations and the functions 
B{x) and C[x ) follow by consistency from the relations 
@. Inserting 03 into 03. and taking into account 
the matter symmetry, Tj? = T®, from the combination 
Rt = R x x it follows that £(x) =constant. Through some 
redefinitions, and without loss of generality, the line ele¬ 
ment m can be rewritten as 

ds 2 = —A(x)dt 2 + } dx 2 + x 2 d0 2 . (18) 

A(x) 

Introducing the standard mass ansatz in three dimen- 
sions, A(x) = —M(x), from the component (g ) of the 
field equations m one finds 


where q is an integration constant identified as the elec¬ 
tric charge associated to a given solution. To formu¬ 
late the matter side of the field equations 03 in a way 
independent of the metric, we note that in terms of 
X = q 2 /{r 2 ip 2 x ) the energy-momentum tensor reads 


T u = ( (‘T’ — 2X ( px)l2x2 0ix2 
\ 62x1 V 


( 14 ) 


M x 




l + ~r — 
y/2 r 


(19) 


The relation between the coordinates x and r(x) also 
follows from Eqs.© as x = fnr(x), which simply im¬ 
plies a shift in the radial coordinate so dx/dr = 1. 
Let us first assume a > 0. Introducing a new variable 
z = y/2r / (3ar g ), we perform the integration of the mass 
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Figure 1. Behaviour of the metric function A (r)//^ for 
f(R) gravity with a Maxwell field (<3 > 0). In this fig¬ 
ure, Mo = r q = 1. From the GR solution ( a = 0, 
dashed curve), the set of solid curves corresponds to grow¬ 
ing values of a = 0.005,0.02,0.05,0.1,0.5. They approach 
the (non-asymptotically flat) GR curve for large r, A(r) ~ 
—Mo — 2r q log[r] + ba/r + 0(a 2 ), with b(Mo, r q ) a constant. 
Note that at r = 0 all solutions with 5 ^ 0 go to zero. For 
a > d c (M = r q = 1) ~ 0.3 the horizon disappears. 


function as M = M 0 + 2 r 2 G(z) where the function G{z) 
satisfies: 


dG l (l + £) 
dz + 


( 20 ) 


Using again Eq. m, we find B{x) = A(x)/fj t and, there¬ 
fore, a full solution to this problem is achieved. For large 
distances, we get fn ~ 1 and dG/dz ~ 1 /z, which upon 
integration yields G(z) ~ log(z), and the standard solu¬ 
tion found in GR, A(x) = —M 0 — 2 r 2 log [ 2 ], is recovered 
0- Note that this asymptotic behaviour is neither flat 
nor (Anti-)de Sitter, which is a generic feature of any such 
solutions in three dimensions in absence of cosmological 
constant term 

The integration of (l20l) is analytical, and can be ex¬ 
pressed as 


G(z) — - + 31og[l + 2 ]^ . (21) 

As one gets close to 2 = 0, deviations with respect to 
the GR solution are found. Going back to the coordinate 
r to compare with the GR case, the function A(r)/f 
can be expanded in series around r = 0 (defining a new 
constant a = 3a/y/2 ^ 0) as 


compared to A GR — > +00 in the GR case) due to the 
factor l//j|, which goes to 0 due to the unboundness 
of the Maxwell field. In terms of horizons, a numerical 
analysis determines that, depending on the values of M 0 
and r q , there will be a single horizon (like in GR) or none 
(see FigHJ). The horizon disappears for a fixed value of 
r q , when a > a c (r q ), where a c (r q ) is a critical value that 
grows with r q . This dramatic change in the behaviour 
of the metric at the center for any a yf 0 translates into 
a worsening of the degree of divergence of the curvature 
invariants. For example, the Kretchsman scalar, K = 
R a , behaves as ~ r q /r 4 in the GR case, but as 
~ 1/r 8 in the f(R) case. This space-time seems to be 
quite pathological. 

Let us consider now the case with a = —a > 0. In 
terms of a variable z = r/(ar q ), we now find G z = 

h )—itt , whose integration yields 


G (z) = ^ ((T~p) + 3lo §[ 2 “ ■ ( 23 ) 

This function ceases to be defined at 2 = 1. Therefore, 
the metric function A(r )/is only defined for r > r q a, 
and behaves around this region as 


A(r)/f 2 R * 
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3 (r — r q a ) 3 


O 


(r - r q a ) 2 


(24) 


which means that for any a / 0 it goes to +00 as 
r —» r q a , which is the same behaviour as the GR solu¬ 
tions. In this case the strength of the central divergence 
is softened, from K = 12 r q /r 4 in the GR case, to the 
dominant term K — 64r^/(d 2 (r — r q a)) 2 . A horizon ex¬ 
ists in all cases, and all solutions approach the GR one for 
r —> 00 (see FigO- These results are quite similar to the 
behaviour of solutions found in four dimensional metric- 
affine f(R) gravity with a nonlinear electromagnetic field 
0 - ' 

Despite the different methodologies pursued in the 
metric-affine and metric approaches for f(R) gravity 
mentioned above, let us mention that four dimensional 
spherically symmetric space-times have been largely in¬ 
vestigated in the latter, including analytical solutions of 
f(R) = R 2 [Hi, Reissner-Nordstrom-like solutions with a 
de Sitter center in f{R) = R + aR 2 with non-linear elec¬ 
trodynamics [36], perfect fluid solutions [37j|, spherically 
symmetric stars and Schwarzschid de Sitter solutions in 
f{R) = R — 'y/R [38, 39|, and constant curvature solu¬ 
tions in several f(R) models flc| . 


[JL + Jl 


2 log[c 


+ --^)r* + 0(r*) 

( 22 ) 

This expansion implies that, for any value of a ^ 0, at 
the center all the solutions approach A(r)//|. —> 0 (as 


III. BORN-INFELD GRAVITY 

We turn now our attention to an extension of GR that 
has attracted a great deal of attention in the last few 
years, dubbed as Born-Infeld gravity [4lj |. with many 
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Figure 2. Behaviour of the metric function A(r)//j| for f(R) 
gravity with a Maxwell field (d = —a > 0). In this fig¬ 
ure, Mo = r q = 1. From the GR solution (d = 0, dashed 
curve), the set of solid curves corresponds to growing val¬ 
ues of d = 0.005,0.02,0.05,0.1, and approach the (non- 
asymptotically flat) GR curve for large r, having a horizon 
always. Note that for any d ^ 0, no solution reaches r = 0, 
but instead they are only defined beyond r > r q a. 

applications in astrophysical and cosmological scenarios 
[42j |. For conceptual and technical simplicity, we can con¬ 
veniently express its action as 

S = -j- [ d 3 x (25) 

where q is the determinant of the new metric q „„ = g^ v + 
eR^, v { r), and e is a small parameter with dimensions of 
length squared, which controls the deviation from GR. 
The meaning of A follows from a series expansion in e as 

S = -4 J d 3 Xy/^g[R - 2A] (26) 

~ 2k? j (—J + 0(e 2 ) 

where the hrst term corresponds to GR with a cosmo- 

,2_i 

logical constant term A = , and the second term 

are quadratic corrections. This just a particular case 
of a family of Lagrangians that can be written under 
the form f(R,R^R ,J ' l/ ). Let us stress that, when for¬ 
mulated in the metric-affine approach, its field equations 
are also second-order and ghost-free, as has been shown 
in four-dimensional settings j43| . It is also worth pointing 
out that a three dimensional model of this f(R, R llv R lil ') 
form has been found a few years ago by Bergshoeff, Hohm 
and Townsend, with the interesting property that its for¬ 
mulation in the standard metric approach is ghost-free 
at the tree level [44] . 

The meaning of the new metric q M! , follows from the 
implementation of the metric-affine approach, so inde¬ 
pendent variations of the action (1251) with respect to the 


metric g^ and the connection yields the system of 
equations 

^Lq» v - A = -K 2 eT» v (27) 

V=9 

V£(>/=^r) = 0 (28) 

so the independent connection is compatible with the 
metric q^ v (but not with the space-time metric g since 
V^(-v/— 99 ^") 7 ^ 0). The Born-Infeld gravity action (l25l) 
can be seen as the comparison of two dimensional volume 
elements between the Riemannian metric q^ and the 
non-Riemannian one AAs shown by Katanaev and 
Volovich [45l] , such a comparison provides a natural mea¬ 
surement of the mass of point defects in crystalline struc¬ 
tures, which is computed according to the difference of 
(two dimensional) volume between a defected crystal (as 
described by g M „) and an idealized perfect crystal with¬ 
out defects (defined by g M „). Since point defects break 
metricity r2.il . the action (12511 in metric-affine approach 
with non-metricity seems to be favoured by the physics 
of crystalline structures. 

From Ea. (|28]) we see that the independent connection 
r* is given by the Cliristoffel symbols of which, 
therefore, carries the same meaning as in the f(R) 
case. Now the relation between g and q^ can also be 
expressed through the algebraic transformation 

Qnu = |t|(T -\ a g av ; <T = 4-^“T a " , (29) 

where we have defined a new object, T = Ifll 1 / 2 ^ -1 , with 
the components of the matrix Cl given by = g a ^qp v . 
In terms of this matrix, the Born-Infeld Lagrangian in 
Ea. (l25l) is simply expressed as fsi = \Cl \ 1 ^ 2 [46| and the 
metric field equations become 

IQI 1 / 2 ^- 1 )^ = A S lt v - en 2 T^ (30) 

so once the energy-momentum tensor is specified, G can 
be fully determined. In addition, this equation implies 
that Q only depends on the matter energy-momentum 
tensor, . 

From the definition eR^( T) = q^ — g M1/ , we raise an 
index with the metric q av and use the above equation to 
write the held equations in terms of q^ as 

i?/( 9 ) = ^[W/ + V], (31) 

where the Born-Infeld gravity Lagrangian is expressed as 

( 32 , 

K Z C 
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In summary, Born-Infeld gravity admits a similar rep¬ 
resentation as f(R) theories in three space-time dimen¬ 
sions. The relation between g^ and v is not conformal, 
but still governed by algebraic transformations only de¬ 
pending on the matter-energy sources, . In addition, 
the same considerations regarding the second-order field 
equations and the recovery of GR in vacuum holds also 
in this case. 

For an electromagnetic field (H4l) , the equation (l30l) can 
be solved by introducing the ansatz 

( 33 ) 


A/Q+ 



where 72x2 and 6 are the identity and zero matrices, 
respectively, and we have defined the objects ft- = 
(A + X) 2 , fl+ = (A — X)(X + X) with X = $£x. The 
field equations (IdTl) in this case take the form 



As in the f(R) case, the source symmetry T t 4 = T x x 
allows to write a metric in the form (frop . and the mass 
ansatz in three space-time dimensions, A(x) = —M(x), 
yields a single independent equation to be solved: 


eM x = x ————. (35) 

From (l29l) and the equations above, it follows that the re¬ 
lation between coordinates in the q and g^ geometries 
is given by x 2 = r 2 fl_, which implies dx/dr = fl + /Q 1 / 2 . 
Plugging this result into (l35l) we obtain 


eM r 


r(il_ — 

nZ 


(36) 


Let us now introduce the length scale e = l 2 and the 
charge scale r 2 = K 2 q 2 /(4:ir) which, together with the 
space-time mass, Mo, characterize the function A(x) = 
— (Mq + G(r)/l 2 ). The function G(r) satisfies 


dG 

dr 



and admits an analytic integration, namely 


(37) 


G(r) = 


1 7 Z 4 r 4 l 2 r 2 

-RiS 1+r2(A2 - 1)+2 — 


log [Z 2 r 2 + Ar 2 ] 


A 


(38) 


Figure 3. Behaviour of metric function A/Q+ for Born-Infeld 
gravity with Maxwell field with A = 1, parameters Mo = r q = 
1, and length scale l 2 > 0. Starting with GR solution ( l 2 = 
0, dashed curve) we increase l 2 = 0.1,0.3,0.5,0.75,1 (solid 
curves), which makes the radius of the core to increase as well. 
A change in the behaviour occurs at l 2 = l 2 c ~ 0.334. Those 
solutions with l 2 < l 2 c always have a horizon, while those with 
Z e > Z 2 c do not. All solutions converge asymptotically to the 
GR case. 


which completely determines the metric function B{r) = 
A /The expansion of this function for r —> oo yields 


A _1-A 2 2 AM 0 +r 2 (21og[r]+log[A]) 

n+ ~ 2Z £ 2 A 2 r A 3 + \r) 

(39) 

which means that we find asymptotically de Sitter so¬ 
lutions if l 2 > 0 and A > 1, or if l 2 < 0 and A < 1, 
while asymptotically Anti-de Sitter solutions are ob¬ 
tained when l 2 > 0 and A < 1, or l 2 < 0 and A > 1. 
When A = 1, the asymptotic behaviour is neither (Anti- 
)de Sitter nor flat, but the standard three dimensional 
GR behaviour, namely, A(r) = — M 0 — 2r 2 log[r]. 

To study the behaviour around the center of the so¬ 
lutions let us first focus on the case Z 2 > 0. We note 
that, given the fact that fl + vanishes at a non-vanishing 
(core) radius r = l e r q / A 1 / 2 , the series expansion around 
it yields the dominant term 


A ~ l e r q (2MpX + r 2 (2A 2 — 1 + 2log[2Z 2 r 2 ])) 
n + ” 8AV2 (r-gfr) 

(40) 

This provides an involved interplay between the mass and 
charge parameters characterizing a solution, Mo and r q , 
the constant characterizing the far-distance behaviour, 
A, and the length scale l 2 . To compare our results with 
those of GR, let us focus on the case A = 1. Given 
the presence in this expansion of the term log[2Z 2 r 2 ], we 
expect a change in the behavior for a sufficiently large 
value of l 2 . Numerical investigation of this question is 
plotted in Figl3] The behaviour at the center depends 
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Figure 4. Behaviour of metric component A/Q.+ for Born- 
Infeld gravity with Maxwell field with A = 1, parameters 
Mq — r q — 1, and length scale l 2 = —l 2 > 0. Start¬ 
ing with GR solution ( l 2 = 0, dashed curve) we increase 
ll = 0.1,0.3,0.5,0.75,1 (solid curves), which makes the ra¬ 
dius of the core to increase as well. All solutions converge 
asymptotically to the GR case. 

on a critical value 1 2 = l 2 c (which depends on the specific 
amounts of mass and charge), in such a way that those 
solutions with l 2 < 1 2 have a single cosmological horizon 
(like in the case of GR), while those satisfying l 2 > l 2 c 
do not, though all cases asymptotically approach the GR 
solution, as already mentioned. 

Let us consider now the behaviour at the center for 
those solutions with l 2 = —l 2 > 0. Now the expansion 
around r = l^q/X 1 ^ 2 becomes 


finite-radius core beyond which they all cease to exist 
grows as well. As in the GR case, there is always a cos¬ 
mological horizon for any value of l 2 (see FigQJ. 

Note that the reason for the existence of this finite-size 
core is due to the presence of the object f2+ in the metric 
function A/Vt + , which vanishes at r = l t r q / A 1 / 2 when 
l 2 > 0, or at r = Urq/X 1 / 2 when l 2 > 0. The latter also 
corresponds to x = 0 due to the change of coordinates 
IPP- Note that, in the equation defining the relation 
between coordinates, x 2 = r 2 f2_, we have extracted a 
square root using the (+) solution. To extend the so¬ 
lution to the region x < 0 in the case l 2 > 0 we could 
consider the (—) solution there, thus replacing x by its 
absolute value in (H51) . Bearing this in mind, the solutions 
could be extended to the region x < 0 using two charts, 
namely, r+ for x > 0 and r_ for x < 0 or, instead, a single 
chart for the coordinate x covering the whole space-time, 
x G] — oo, +oo[. However, the matching of these geome¬ 
tries at x = 0 is not smooth since dr/dx\ x —o ^ 0 [see 
Ea. (H51) ]. Indeed, at x = 0, using Mathematica we find 
the Kretchsmann scalar, K = R a 1 is dominated 

by a curvature divergence of strength 

K ~ log [x] 2 /x 2 , (44) 

which is milder than in the GR case (for large distances 
we get K ~ 12r 2 /r 4 , the standard GR behaviour). De¬ 
spite this unpleasant behaviour, the surface x = 0 have 
some interesting properties. For example, an explicit 
computation of the electric charge flowing through the 
surface x = 0 can be done through the flux integral 


A l e r q (2M 0 X + r 2 (l - 2X 2 + 21og[r 
n + 8A 7 / 2 (r - £5*) 

+ 2f e r g log[2 lr q \W] + 

8A 7 /2 (r - fa) 



D) 


(41) 


Let us focus again on the case A = 1. The line element 
for g^ v can now be conveniently expressed as 


$ = / *F = 2nq, (45) 

Js 1 

where *F is Hodge dual, and the integration takes place 
over a boundary defined by close hypersphere S 1 , where 
2-7T is the volume of the unit S 1 sphere. An important 
property of this flux is that its rate of growth with the 
core “area”, S = 2nr , defines a natural constant, i.e., the 
density of the flux 


ds 2 


A 2 

-FT dt +^~ 
r2_i_ A 



+ r 2 (x)d0 2 , 


(42) 


where the radial coordinate r(x) is explicitly written as 


r( x) 


x+^/x 2 + 4r 2 l 2 
2 


(43) 


The behaviour close to the center, for any values of 
the parameters M 0 and r q , and length scale l 2 is qual¬ 
itatively the same. For growing values of l 2 each curve 
smoothly deviates from their GR counterparts, and the 


JL = (46) 

27 rr l e 

is independent of the particular amounts of mass and 
charge, depending only on the fundamental constants G, 
c (through Newton’s gravitational constant in three di¬ 
mensions, k 2 ) and the length scale l e . When l 2 —► 0, 
this finite-structure disappears and the standard point- 
like singularity of GR is recovered. However, the non¬ 
smoothness of the change of coordinates dr/dx at x = 0 
seems to prevent a wormhole-type extension, as opposed 
to the case of four dimensional Born-Infeld gravity [43j . 

Let us now discuss briefly the case of solutions with 
A / 1, according to the asymptotic behaviour given by 











Ea. (l39l) . We have investigated numerically this question 
and the following structures are found: 

When ll > 0 and A > 1 we find the same num¬ 
ber and type of horizons as in the GR case, but now 
with an asymptotically de Sitter behaviour, instead of 
~ — 2r^log[r]. For > 0 and A < 1 the behaviour is 
asymptotically Anti-de Sitter and we find that, as in the 
cases with A = 1, a critical value 1% (M 0 , r q , A) determines 
a sudden transition. In this sense, those solutions with 
11 < correspond to naked singularities, while those 
with 11 > are black holes with a single non-degenerate 
event horizon in all cases. 

When ll > 0 we find asymptotically de Sitter solu¬ 
tions corresponding to A < 1, which have a cosmologi¬ 
cal horizon in all cases. However, for A > 1 we find an 
asymptotically Anti-de Sitter behaviour and three classes 
of structures, depending on a critical value (Mo, r q , A). 
In this sense, solutions with 11 > /jy correspond to black 
holes with two horizons (inner and event), which merge 
into a extreme, degenerate horizon (extreme black holes) 
when = 1% , and solutions with l\ < lj c correspond to 
naked singularities. 

These behaviours resemble what it is usually find in 
other (Anti-)de Sitter backgrounds, like in the context of 
Gauss-Bonnet gravity |47|. 

IV. CONCLUSIONS 

In this work we have considered two families of mod¬ 
ified gravity theories sourced by electromagnetic fields. 
Using results from crystalline structures in solid state 
physics, where the presence of point defects on their mi¬ 
crostructure requires a metric-affine framework for their 
proper description at the macroscopic level, we have ob¬ 
tained the equations for these theories from independent 
variations with respect to metric and connection. Moti¬ 
vated by the existence of some ordered structures with 
defects which can be effectively treated as bi-dimensional 
systems (like graphene), we have considered metric-affine 
gravities with non-metricity in a three dimensional space- 
time. The resulting formalism has several advantages 
with respect to the standard metric formalism of mod¬ 


ified gravity, like second-order equations in all cases or 
the recovery of GR in vacuum. 

Explicit solutions to an f(R) model were found and 
characterized. While for large distances they boil down 
to their GR counterparts, close to the center large devia¬ 
tions are found. The case of f(R) with a > 0 seems to be 
quite pathological, as the metric at the center of the so¬ 
lutions vanishes for any finite value of a , with a worse de¬ 
gree of divergence of the curvature invariants there. This 
has led us to consider a crystal-motivated action, as given 
by Born-Infeld gravity, where a finite-size structure arises 
as long as l\ is non-zero, replacing the standard point¬ 
like singularity of GR. When A ^ 1, and depending on 
the choice of parameters, standard horizon structures of 
black holes may arise. Our results support the interpreta¬ 
tion of 11 as a natural scale around a point defect where 
the geometry may largely deviates from the GR coun¬ 
terpart. Further investigation on these geometries could 
avoid the shortcoming of the typical space-time singular¬ 
ities of the Riemannian structure which, however, seems 
to be in conflict with the fact that well defined geometries 
exist in the case of crystalline structures with defects. To 
get further into this question, one should consider the 
geodesic behaviour in these scenarios with more realistic 
energy-matter sources, a topic which has received little 
attention in the literature so far. In summary, the metric- 
affine formulation seems to be a more appropriate for¬ 
malism to deal with defected crystals, which encourages 
further research on the field to identify potential experi¬ 
mental signatures resulting from non-metricity. Progress 
in this sense is currently underway. 
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